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Abstract. High-frequency quasi-periodic oscillations detected in the light curves of black hole can- 
I didates can, according to one model, be identified with hydrodynamic oscillations of the accretion 

Odisc. We describe a non-linear coupling mechanism, suggested by Kato, through which inertial 
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waves trapped in the inner regions of accretion discs around black holes are excited. Global warping 
and/or eccentricity of the disc have a fundamental role in this coupling: they combine with trapped 
modes, generating negative energy waves, that are damped as they approach the inner edge of the 
disc or their corotation resonance. As a result of this damping, inertial oscillations are amplified. 



We calculate the resulting eigenfunctions and their growth rates. 
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1. INTRODUCTION 



The idea of using wave coupling as an excitation mechanism for oscillations in de- 
formed discs goes back to Goodman [1], who studied the excitation of waves in tidally 
deformed discs via the parametric instability. Papaloizou & Terquem [2] also mention 
the possibility of parametric generation of inertial waves in warped discs. More recently, 
an analytical study of the excitation of waves in warped and eccentric discs around non- 
rotating black holes was carried out by Kato [Hi 3]. His estimates for the growth rates of 
several modes in the disc are interesting, but many aspects of his calculations are uncer- 
tain, and the origin and dynamic properties of the global warping or eccentricity are not 
discussed. 

In the work reported here, Kato's ideas are developed and generalised for discs around 
rotating black holes, and detailed numerical calculations of trapped inertial modes and 
their growth rates are performed. A dynamical treatment of global deformations is 
included; their origin is discussed in Ogilvie & Ferreira (this volume). 

In §2 we discuss the trapping of inertial waves, in §3 we describe the coupling 
mechanism through which these waves are excited, in §4 we discuss the variation of 
the inertial mode growth rate with several parameters, and we present our conclusions 
in §5. 

2. TRAPPED OSCILLATIONS 



Trapping of oscillations [|5|] is the result of the non-monotonic variation of the epicyclic 
frequency K in relativistic discs. As the radius decreases, K increases, reaching a maxi- 



mum, and then decreases rapidly, going to zero at the marginally stable orbit (regarded 
as the inner edge of the disc). Analysis of the hydrodynamic equations in a simple disc 
model [6] shows that inertial waves can be trapped below this maximum [Q] [see also 



Here we consider an isothermal disc with a ratio of specific heats 7 = 1 . We ignore 
viscosity and magnetic fields, and relativistic effects are included by using relativistic 
expressions for the angular velocity Q., vertical frequency £l z , and K ft 1 Oil . Perturbations 
about the equilibrium disc can be described by 

(u' r , wj, , ti) = Re I [ur(r) , uq (r) ,h(r)] He„ (J^j exp(im^) - icot) j , (1) 

u' z = Re jw z (r)He„_i exp(im0 — icw)j , (2) 

where (u' r , uL,u' z ) are the components of the velocity perturbation, h' is the enthalpy per- 
turbation, m and n = 0, 1,2, ... are the azimuthal and vertical wave numbers, respectively, 
CO is the oscillation frequency, and H = c s / Q. z is the vertical scale-height of the disc, with 
c s being the constant sound speed; (r,0,z) are cylindrical polar coordinates. The radial 
structure of linear oscillations is then described approximately by the following set of 
ordinary differential equations in r 
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where (b — CO — m£l is the Doppler- shifted wave frequency, which is zero at the corota- 
tion resonance. 

The dispersion relation for waves with local radial wavenumber k, can be obtained 
from these equations by assuming the radial wavelength of oscillations to be much 
smaller than both the azimuthal wavelength and characteristic radial scale in the equi- 
librium state; it reads 

2 _ (d) 2 -K 2 )(d) 2 -^ 2 ) 

From this relation we can determine the different types of waves existing in a disc, 
as well as their propagation regions. If n 7^ 0, acoustic (p) modes can propagate where 
a 2 > max(?c 2 ,rc£2 2 ) = n£2 2 , and inertial (r) modes with cb 2 < min(fc 2 ,n^ 2 ) = k 2 can 
exist. If n = 0, waves can propagate where cb 2 > fc 2 ; this is the inertial-acoustic mode. 

Here we focus on the simplest possible inertial mode (m = 0, n = 1 and simple 
radial structure), with frequency CO ~ max(fc), which is confined within a small region 
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FIGURE 1. Variation with radius of the real part of the radial velocity of the simplest possible r mode 
(l,m,n) = (0,0,1), for c s /c — 0.01, a = 0.5. The maximum is achieved at the radius where K is also 
maximised. The crosses and asterisks indicate the Lindblad (ft) 2 = K 2 ) and vertical (co 2 — Q. 2 ) resonances, 
respectively. 



below the maximum of the epicyclic frequency, making it more likely to survive in the 
presence of turbulent viscosity. An axisymmetric mode is likely to be the most relevant 
observationally, as it results in no cancellations in the net luminosity variations of the 
disc. The fact that the frequency of this mode can be identified with max(K") is important 
if one is interested in studying the properties of the black hole, since the maximum of 
the epicyclic frequency depends only on the mass M and spin a of the compact object 
tfH. 

To find the radial structure of axisymmetric modes with n = 1, we solve equations ©- 
© numerically QVU as a generalised eigenvalue problem, subject to appropriate boundary 
conditions. (For more details on the numerical method and boundary conditions used see 
I1120-) We find the simplest possible (lowest radial order, / = 0) trapped mode to have an 
approximately Gaussian structure in r, centred at the maximum of K, in accordance with 
UM (see Fig. [D. 



3. EXCITATION MECHANISM FOR TRAPPED R MODES 

We are interested in studying the growth of trapped inertial oscillations in deformed 
discs. Deformations like warping or eccentricity are typically identified with low- 
frequency modes with m = 1, as these oscillations are global, long-lived, and vary on 
length- scales comparable to the characteristic disc radius. The warp is identified with 



an (m,n) = (1,1) mode jl4\. while eccentric discs correspond to m = l,n = modes. 
Around rotating black holes, the stationary warp and eccentricity have an oscillatory ra- 
dial structure [15] (see also Ogilvie & Ferreira, this volume). The excitation mechanism 
for r modes relies on their interaction with these global deformations, which gives rise 
to an intermediate mode that can then couple with the warp or eccentricity to feed back 
on the original inertial oscillations. 

The basic hydrodynamic equations have an intrinsic non-linearity through the term 
u • Vu, which provides couplings between the different linear modes of the system. 



Here we focus on the couplings between the global deformations, and the inertial and 
intermediate modes that result in excitation of the simplest possible trapped r mode. 

Before describing our numerical calculations on the excitation mechanism, we should 
analyse the energy exchanges between the different interacting modes, to make sure the 
r mode gains energy in this coupling. We consider the warp or eccentricity to have zero 
frequency so that it propagates everywhere in the disc (according to the dispersion re- 
lation), and therefore its energy is essentially zero. Consequently, the energy exchanges 
happen between the r and intermediate modes and the disc. Since the first is axisym- 
metric, it has positive energy, while the intermediate mode has negative energy, as it 
propagates inside its corotation radius. In the coupling between the global deformation 
and the r mode, the latter is amplified in the process of generating the negative energy in- 
termediate wave. Sustained growth of the inertial mode requires this intermediate wave 
to dissipate, drawing positive energy from the disc rotation, so that its negative energy is 
continually replenished by the r mode. To achieve this dissipation, we choose to damp 
the intermediate wave locally at a rate /3Q, where /3 is a dimensionless parameter, i.e. 
we include a friction term in the equations for the intermediate mode. This does not 
necessarily mean that the mode is dissipated into the disc due to viscosity or friction. It 
is also a way of representing mathematically the fact that the intermediate mode might 
be absorbed at its corotation radius or at the marginally stable orbit. 

For wave interaction to occur, the modes need to propagate in the same region of 
the disc, and the parameters (ft),m,n) need to obey certain coupling rules. For the 
interaction of the r mode ((0,0, 1) with the warp (0, 1, 1), two intermediate modes are 
possible: (ft), 1,0) and (ft), 1, 2). The interaction with eccentricity (0, 1,0) results in only 
one intermediate mode (ft), 1, 1). For each one of the three possible interactions through 
which the inertial mode might be excited, we solve a system of equations for the r and 
intermediate modes, coupled by the warp or eccentricity J12Q. Even though the coupling 
terms arise from the non-linearities in the basic equations, the system is linear in the 
unknowns (the enthalpy and the components of the velocity of the r and intermediate 
modes). Expressions for warp or eccentricity enthalpy and velocity perturbations are 



Each of the systems of coupled equations is solved numerically as an eigenvalue 
problem, subject to certain boundary conditions. We aim to find solutions for which the 
radial velocity of the r mode resembles the function shown in Fig. [H i.e. we choose the 
coupling terms to be small enough that the eigenfunctions are not significantly altered, 
and the eigenvalues (— ift)) change slightly by acquiring a small positive real part. The 
growth rate of the inertial mode is then given by the (positive) imaginary part of its 
frequency. 

In Figs [2] (a) and (b) we show the variation of the growth rate with the warp amplitude 
at the inner radius, and with the sound speed in the disc, respectively, for the interaction 
where the intermediate mode has n = 0. Results are similar for the interaction with the 
n = 2 mode, as well as for the interaction in the eccentric disc, in which case the warp 
amplitude at the inner radius is replaced by the eccentricity at the same location. In Figs 
|2] (c) and (d) we represent the variation of the growth rate with the spin of the black 
hole for the interaction of the r mode with the warp and n = mode (results for the 
n = 2 mode are similar), and for the interaction in the eccentric disc. These results are 
discussed in the following section. 
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FIGURE 2. Panels (a) and (b) show the variation of the growth rate of the simplest trapped r mode, 
with warp amplitude at the inner radius, and sound speed in the disc, respectively, for the interaction with 
the n = intermediate mode. Panels (c) and (d) show the variation of the growth rate with the spin of the 
black hole, for the interaction with the n — intermediate mode in a warped disc, and the interaction with 
the n — 1 intermediate mode in an eccentric disc, respectively. 



4. DISCUSSION 

In Fig. [2] (a) we can see that the growth rate of the inertial mode varies with the square 
of the disc inclination at the inner radius. This is expected since the interaction relies on 
the use of the warp twice, and we can therefore argue that the growth rate is proportional 
to |dW / dr| 2 , where W(r) is the inclination of the disc at radius r (and dVK/dr represents 
the actual warp). In turn, |dW/dr| 2 is proportional to the square of the inclination at the 
inner radius divided by the local wavelength of the warp (Xw) squared, which, according 
to the dispersion relation for (co,m,n) = (0, 1, 1), is proportional to c 2 . Therefore, for a 
constant sound speed, we expect the growth rate to be proportional to the square of the 
inclination at the inner radius (Fig. [2] (a)). Also, if the latter is fixed, the growth rate is 
proportional to the inverse of the square of the sound speed, as we can see in Fig.[2](b). 

The variation of the growth rate with the spin of the black hole for the interaction 
in the warped disc (Fig.[2](c)) shows that there is no mode excitation if the black hole 
is non-rotating. Although this result is in apparent contradiction with Kato's work 
we can argue that it is in fact expected. If the black hole is not rotating, the metric 
is spherically symmetric, and the warp is only a rigid tilt, which means that nothing 



changes in the disc except its inclination angle. If there are no relevant changes in the 
disc, there should be no relevant changes in the oscillations, i.e., a rigid tilt warp should 
not work as an excitation mechanism for waves in a disc. However, even in the case of 
a non-rotating black hole, it is still possible to have growth of oscillations, if the disc 
around it is eccentric (Fig.[2](d)). 

5. CONCLUSION 

We have described an excitation mechanism for inertial modes trapped in the inner 
region of discs around black holes. The mechanism relies on a non-linear coupling 
between these waves and eccentricity or warping of the disc. Excitation of the inertial 
modes is possible if the intermediate mode resulting from the interaction deposits its 
negative energy into the disc, by being dissipated or absorbed at its corotation radius 
or the marginally stable orbit. Rotational kinetic energy is then removed from the disc 
and becomes available for the r mode to grow. Depending on the disc properties (sound 
speed and amplitude of the warp or eccentricity at the inner radius) and on the spin of 
the black hole, reasonable values of the growth rate can be achieved for a deformation 
of modest amplitude. 
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